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Abstract 

One-prong r decays into final states involving kaons are studied with about 
161k r+T~ events collected by the ALEPH detector from 1991 to 1995. Charged 
kaons are identified by dE/dx measurement, while i^^^'s are detected through their 
interaction in calorimeters. Branching ratios are measured for the inclusive mode, 
B{t~ K^Xvr) = (1.52 ± 0.04 ± 0.04)%, where X can be any system of neutral 
particles, and for the exclusive modes 
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(6.96 ±0.25 ±0.14) x 10^^ 
(4.44 ±0.26 ±0.24) x lO'^ 
(0.56 ±0.20 ±0.15) X 10-3 
(0.37 ±0.21 ±0.11) X 10-3 
(1.62 ±0.21 ±0.11) X 10-3 
(1.43 ±0.25 ±0.21) X 10-3 
(9.28 ± 0.45 ± 0.34) x 10-3 
(3.47 ±0.53 ±0.37) x 10-3 



Upper limits for B(t K^t:' and B{t K K^Ti'^Tx'^Vr) are also ob- 

tained. Mass spectra in the final states are addressed in order to study the relevant 
dynamics. 
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1 Introduction 



The study of r decays involving kaons is necessary in order to understand the strange 
hadronic r sector as such, but also to provide a better control over nonstrange decay 
modes. Both issues play an important role in the low energy QCD analyses accessible 
through hadronic r decays. 

Results on charged kaons in three-prong decays [ |I| and on production [ |^ 
have already been published. In this paper, one-prong decays are investigated. Taking 
advantage of their simpler topology, the analysis combines the dE/dx measurement of the 
charged track and the reconstruction of -ft'^'s through their hadronic showers. 

The first analysis of one-prong final states with a charged kaon was presented using 
38k selected r+r" events collected in 1991 and 1992 [ and an update was given adding 
the 1993 data for a total of 64k selected r+r^ events [ ^, enabling the study of the 
inclusive kaon mode and the three exclusive channels^ K~Vr-, K'tt'^Ut-, and K^tt^tt^Ut-. 
Analyses with detection were published with the same data sets [ |^. In this 
paper, final results are given using the full LEPl statistics corresponding to 161k selected 
r+r~ events. Several improvements are implemented, increasing selection efficiencies and 
reducing systematic uncertainties. This is achieved through a better use of all the available 
calorimetric information, imbedded into a more efficient method for measuring the K^^ 
signal. Because of the larger statistics, a better understanding of the systematic effects 
is also acquired. The present results therefore supersede those from previous analyses on 
one-prong kaon modes [ |, |[ |]. 

After a short description of the ALEPH detector and the event selection, the dE/dx 
calibration using muons is presented in detail. Then, the calorimetric method for 
detection, taking advantage of the granularity of the calorimeters and the special r decay 
kinematics, is described. Branching ratios for inclusive charged kaon production and 
exclusive channels involving charged and neutral kaons are derived through statistical 
methods, both for dE/dx and calorimetry. A detailed account of the systematic 
uncertainties is given next, showing that in all cases the measurements are still statistically 
limited. Finally, the underlying dynamics in these decays is addressed through the study 
of invariant mass spectra for exclusive decay samples. 

2 The ALEPH detector 

A detailed description of the ALEPH detector can be found elsewhere [ ^ . Here only the 
most relevant features for this analysis are described. 

The momenta of the charged tracks are measured by means of a central tracking system 
consisting of a silicon microstrip vertex detector, an inner tracking chamber (ITC), and 
a time projection chamber (TPC). All of these three tracking devices are immersed in a 
1.5 T axial magnetic field, and achieve a transverse momentum resolution for the charged 
track of apjpt = 6 x lO'Vt® 0.005, with pt in GeV/c. 

Two subdetectors, the electromagnetic calorimeter (ECAL) and the hadron 
calorimeter (HCAL), produce the energy measurement for electrons and hadrons. The 

^Throughout this paper, charge conjugate states are imphed. 
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ECAL consists of 45 layers of sandwiched lead sheets and proportional wire chambers, 
representing 22 radiation lengths, and is used to measure the energy deposit of electrons 
and photons. Its fine granularity is essential to the reconstruction of vr'^'s, and the 
longitudinal segmentation with three stacks measures the electromagnetic profile of the 
detected shower. The energy resolution in the ECAL is measured to be (Te/E = 

.18 / \J E (GeV) + 0.09. Hadronic energy is measured by the HCAL, which is composed of 
a 1.2 m thick iron stack serving as the yoke for magnetic flux return. Streamer tubes are 
inserted into the gaps of the iron stacks, giving 23 energy samples. The energy resolution 
of the HCAL for pions at normal incidence is <7e/E = 0.85/ ^E{GeV). Outside the 
HCAL, two layers of streamer tubes are used to complete the muon detection. 

In the data sample used for this analysis, all components of the detector are required 
to be operational. Furthermore, one of the following trigger conditions is required: a 
minimum ECAL energy of 6 GeV, or a track segment in the ITC pointing to an energy 
deposit in the ECAL of at least 1.2 GeV, or a track segment matching the signal of a 
penetrating particle in the HCAL. By comparing redundant and independent triggers 
involving the tracking detectors and the main calorimeters, the trigger efficiency for r 
pair events is estimated to be better than 99.99% within the selection criteria used in this 
analysis. 

The measurement of the ionization dE/dx deposited by the charged particles when 
traversing the TPC volume is the only tool in ALEPH for the identification of charged 
hadrons. Ideally, the dE/dx ionization can be measured with up to 338 samples. However, 
in some case a charged track may not have such a number of samples because of the 
following factors: the cracks between the TPC sectors, the geometrically available track 
length and the applied truncation for avoiding the large fiuctuations caused by 6 rays. 
A minimum of 40 samples is required to ensure a good particle identification, rejecting 
about 0.1% of the events. 

3 Event selection 

A sample of about 161k events is obtained using a general r-pair selection procedure [ 
H, 1^, from about 202k produced r pairs. The one-prong decays are selected from this 
sample. 

Each event is first divided into two hemispheres separated by the plane perpendicular 
to the thrust axis of the event. A one-prong hemisphere is required to contain exactly one 
good charged track, satisfying a minimum of four associated hits in the TPC, an impact 
parameter in the plane transverse to the beam axis |cio| < 2 cm and the distance from the 
interaction point along the beam axis \zq\ < 10 cm. Also, the track momentum is required 
to be larger than 100 MeV/c. To reduce the number of fake K^^s produced by additional 
tracks, each hemisphere candidate is not allowed to contain any unselected charged track 
associated with at least three hits either in the TPC or the ITC, |(io| < 20 cm, \zo\ < 40 cm 
and a momentum more than 2 GeV/c. Finally, to identify electrons, muons and hadrons 
in the one-prong r candidates, a particle identification [ 0] is performed, except for the 
tracks with momenta below 2 GeV/c (where only electrons are identified through their 
dE/dx) or going through a crack between ECAL modules. All these requirements select 
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a total sample of 143 319 one-prong hadronic r decays. Hemispheres containing a track 
with momentum less than 2 GeV/c and not identified as an electron by dE/dx are kept 
for the search performed in Section 7. 

Then, all selected one-prong candidates are classified according to the number of 
associated vr^'s, which are identified using a photon and 7r° reconstruction algorithm [ 
In this analysis, a vr*^ is defined as follows. A resolved vr^ is formed by a pair of 
photons with an invariant mass consistent with that of a vr''. Photons from high energy 
TT^'s may give rise to electromagnetic showers so close to each other that they cannot 
be recognized as two separate objects, leading to an unresolved ti^ which is identified 
by the energy- weighted moment technique [ The remaining photons inside a cone 
of 30° around the thrust axis are called residual photons. Under these definitions, one- 
prong r decays are separated into h~ {h = n or K), h^n^, /i~7r°7r°, /i~7r°7r°7r° and the 
rest (containing one or several residual photons), corresponding to 47543, 54171, 10 325, 
959 and 30 321 decays, respectively. 

Finally, to evaluate the efficiencies and the background contamination, 1.0 million 
r+r~ Monte Carlo events and 5.8 million qq are produced by the standard KORALZ [ 
^ and the JETSET 7.4 [ ^ generators, respectively, with a complete ALEPH detector 
simulation. 



The ionization energy deposit R is measured in the TPC The expected Ri and its 
standard deviation cxj for a given particle type hypothesis i are initially calculated using 
a parametrization of the Bethe-Bloch formula fitted to hadronic events [ ^ , in which the 
environment is rather different from that of the one-prong r decay where no track overlap 
occurs. In one-prong hadronic r decays, the environment is similar to that of the muons 
from Z or r~ — > ^'v^v^ which provide a pure sample of isolated muon tracks 

covering the full momentum range and characterized by the same angular distribution in 
the detector as the hadron sample under study. Therefore, a special and accurate dE/dx 
calibration can be performed in order to achieve a more precise measurement on the small 
kaon fraction in r decays. 

The main philosophy of dE/dx calibration follows Ref. [ ^. Some modifications are 
made, taking into account slight variations in each period of data taking and allowing for 
possible small differences between muons and hadrons, due to interactions for the latter. 
The calibration is based on a sample of muons selected from Z — > and 
— s> ^^v^Vri corresponding to 278 079 and 49 620 tracks respectively. Distributions 
of the variable = '~ are studied as functions of (i) the period of data acquisition; 
(ii) the polar angle of the track; (iii) the number of dE/dx samples; and (iv) the R value. 
The calibration gives the new values for R^ and according to 
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Figure 1: distribution for the muon sample (top). The fit is indicated by the solid 
curve, and the expected residual kaons and pions are shown by dashed and dotted curves, 
respectively. The residuals between the data and the fit are shown in the bottom plot. 

where x"^'^ and a^ou represent the mean and the standard deviation of the initial x^ 
distribution. The cahbration is performed in an iterative way correcting most of the 
correlations among the above four factors. The distribution of the calibrated x^ for the 
muon sample is parametrized by an initially energy-independent four-Gaussian function 
F^{x) showing a slight asymmetry. There is still a tiny contamination from residual kaons 
and pions in the selected muon sample from r's because of either hadron misidentification 
as muon or decay in flight. These backgrounds are estimated directly from the Monte 
Carlo simulation, and the corresponding fractions are fixed in the fit. The final result 
of the parametrization to the muon sample is shown in Fig. |1|, indicating an excellent fit 
with x^/ndf = 42.9/46. 

The energy independence of F^{x) is not well satisfied at low momentum, as 
demonstrated with a large statistics sample of 77 — > events. It is observed that, 

in the region below 7 GeV/c, the function becomes more symmetric. This small effect is 
taken into account in the final parametrization of the resolution function. 

Finally, the Xh{h) distribution is treated to be of the same functional form as used for 
muons, namely 

Fhixh) = F^{x^). (3) 
With this function, the expected dE/dx response can be generated for any particle type. 
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Figure 2: Separation from pions in standard deviations for different particle types as a 
function of momentum. 

resulting in a distribution from kaons in terms of of the form 

Gk{x^) = — Fk . (4) 

The final particle separations in standard deviations are illustrated as a function of 
momentum in Fig. 

5 Calorimetric calibration 

Since the component is extracted from the calorimetric response by means of a 
statistical technique, a calibration of the measurements in both the ECAL and the HCAL 
is necessary. The details of the calorimetric calibration are described below. 

5.1 Definition of the discriminating variables Scj) and 5e 

At LEP energies, the produced r leptons are strongly Lorentz-boosted with only a few 
degrees opening angle between any two particles in the final state. However, the bending 
due to the magnetic field provides an opportunity for separating the calorimetric clusters 
associated to the and the charged hadron in the HCAL, by making use of the transverse 
angular offset 

i I (pbarycentre 4^ track impact I, (5) 

where the positive sign is taken if the azimuthal angle (pbarycentre, determined by the 
barycentre of energy measured in the HCAL, is inside the bending trajectory of the 
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charged hadron, and the negative sign is taken otherwise; (ptrack impact is defined as the 
azimuthal angle for the position of the charged track at the inner radius of the HCAL. The 
minimum transverse momentum needed for a charged particle to reach the HCAL is about 
0.7 GeV/c provided that no energy is lost in the ECAL. K^^s from r decay are emitted 
at small angle with respect to the charged hadrons, leading to a characteristic negative 
value of 6(f) in most of the cases. Throughout this paper 6(j) is expressed in degrees. 

In addition, the signal can be enhanced by exploiting the hadronic energy excess 
in the calorimeters, measured by the quantity 



where E is the sum of the energies measured by both the ECAL and the HCAL in the 
studied hemisphere, E^^o is the total 7r° energy measured by the ECAL in hir^lir^) modes, 
Ph is the charged track momentum, the denominator is close to the actual resolution of 
the calorimeters for hadrons, and C = 0.9 (GeV/c)2. Since -R'^'s from r decay have a 
momentum in excess of 3.5 GeV/c, a positive value of 6e is expected. This procedure 
differs from the previous analyses [ ^ which used the energy deposit in the HCAL 
only. In fact more than 50% of K^^s interact inside the ECAL and leave part of their 
energy mostly in the last two stacks, with longitudinal and transverse profiles different 
from those of a photon shower. Since the simulation of the relative energy deposits in the 
ECAL and the HCAL is rather delicate, it is more reasonable to combine the two pieces 
of information, in order to reduce systematic uncertainties. 

In the following, cuts on 6(j) and Se are used separately in order to select samples for 
calibration. 

5.2 Calibration for the energy measurement 

Decays with 6(j) > 0, consisting of single charged hadrons with a very small contamination 
from and unreconstructed tt*^, are used for the energy calibration. The total hadronic 
energy measured by the calorimeters is summed with a 1.3 weight factor for the ECAL, 
taking into account the ratio of the electromagnetic calorimeter responses to electrons and 
pions [ 1^. The calibration provides corrections as a function of the data-taking period 
and the polar angle, the latter corresponding to varying effective interaction length and 
geometry of the calorimeters. The same procedure is performed on the Monte Carlo 
sample so that its 6e distribution matches data. 

The calibration procedures are applied separately for the h, hir^ and hir^n^ samples. 
Figure ^(a) shows a good agreement between data and Monte Carlo simulation in the 
calibrated 5e distributions for all the modes. The above calibration is applied to all 
events, assuming that events with a negative 50 follow the same variation as those with 
a positive 50. 

5.3 Calibration for the barycentre offset Scj) in the HCAL 

At this stage, a non-K^ sample is selected by requiring hemispheres with a negative 
sign oi 5e- The 50 distribution is checked and found to be consistent within different 
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Figure 3: (a) The calibrated 6e distributions for hemispheres with a positive 6(j). (b) The 
calibrated 5(j) distributions for hemispheres with a negative 6e- Data and Monte Carlo 
simulation are shown in dots with error bars and in histograms, respectively. 

periods of data acquisition. For all hadronic samples, the Monte Carlo simulation is in 
good agreement with data for the 6(j) mean, but produces a somewhat larger width. This 
effect is proportional to the energy deposited in the HCAL, and is used as a correction to 
the width of the 6(f) distribution as a function of energy in the Monte Carlo simulation. 
Figure ^(b) shows that for the calibrated 6(f) distribution in the h mode, Monte Carlo 
simulation is observed to agree with data. The same calibration is applied to the events 
with a positive 6e- 

6 Channels involving charged kaons 

In this section, the inclusive one-prong r decay is first investigated to finalize the dE/dx 
calibration, and a measurement of the corresponding kaon fractions is presented. The 
dE/dx measurement is also exploited to study the exclusive one-prong r decays, which are 
classified into -ft^^-reduced and ii'^-enriched samples. Branching ratios are measured for all 
the exclusive modes R-, R-n^, R-tt^tt°, R-ji^ti^ti^, R-R^, R-R^'^ and R-R^ti^ti^. 

6.1 Charged kaon fraction in the inclusive mode 

To be insensitive to the a priori unknown dynamics, the fit to the x,r distribution is 
performed in momentum slices. Small contaminations from residual muons and electrons 
are included in the fit, with fractions estimated by the simulation. The R/tt separation, 
varying with momentum as shown in Fig. ^ is determined by Monte Carlo and kept fixed 
in the fit. The x^(7r) parameters and the kaon fraction are free to vary. Table |l| lists the 
fit results in the different momentum slices. Figures ^ and ^ show the sum of the fits to 
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p (GeV/c) 


— / \ 

X7r(7r) 




fK{%) 


Nk 




2 / AC 
X /ndf 


2.0 — 3.5 


—0.03 lb 0.01 


1.00 ± 0.01 


0.11 ± 0.12(0.11) 


10 ± 10 


6 ± 2 


A f\ A 1 A ^ 

49.4/47 


3.5 - 7.5 


-0.03 ± 0.01 


1.00 ± 0.01 


1.30 ± 0.11(0.10) 


320 ± 25 


6 ± 2 


69.5/57 


7.5 - 11.5 


-0.02 it 0.01 


0.99 ± 0.01 


2.82 ± 0.15(0.14) 


625 ± 31 


4 ± 2 


73.6/53 


11.5 — 15.5 


—0.02 ± 0.01 


0.99 ± 0.01 


3.92 ± 0.23(0.20) 


680 ± 34 


8 ± 2 


49.3/50 




_n 4- n ni 

— u.uo HI u.u± 


n OS 4- n ni 




uu 1 m oo 


1+1 

1 HI L 


QC //I Q 
OO.U / ^C7 


19.5 - 23.5 


-0.01 ±0.01 


0.98 ±0.01 


4.83 ± 0.44(0.32) 


508 ± 33 


7±2 


38.6/42 


23.5 - 27.5 


-0.03 ±0.01 


1.01 ±0.01 


3.48 ±0.50(0.35) 


303 ± 30 


5±2 


63.7/45 


27.5 - 31.5 


-0.02 ±0.02 


0.99 ±0.01 


4.17 ±0.58(0.39) 


303 ± 29 


5±2 


46.3/39 


31.5-35.5 


-0.02 ± 0.02 


0.98 ± 0.02 


4.30 ± 0.69(0.46) 


247 ± 27 


7±2 


42.5/38 


35.5 - 39.5 


-0.02 ± 0.02 


0.98 ± 0.02 


4.35 ±0.99(0.60) 


199 ± 28 


5±2 


38.5/36 


> 39.5 


-0.02 ± 0.02 


0.98 ± 0.02 


5.04 ± 1.21(0.70) 


204 ± 28 


16 ±3 


53.4/39 



Table 1: Charged kaon momentum range, the fitted XT^{7r) parameters and kaon fractions, 
kaon yields, estimated non-T background together with the test of goodness of fit. The 
errors in parentheses are those without the contribution from the uncertainties on the 
^■k{t^) parameters. 



for all momentum slices and the kaon yield as a function of momentum for both data 
and simulation. 

Due to the large statistics, the inclusive fits yield the final calibration constants 
(relative to the muon calibration) as a function of momentum which are used in the later 
fits to the distributions for all the exclusive modes. For the whole momentum range, 
the obtained values = — 0.021±0.004 and = 0.992±0.003 are significantly different 
from those (with similar errors) from the t ^ ji sample which covers the same momentum 
range. These shifts are consistent with those found using a 99.5% pure pion data sample 
in the h^Ti^ mode by requiring the invariant mass M(i^~7r°) > 1.1 GeV/c^, in order to 
veto the K*{m2)- background. The values = -0.02 ± 0.01 and = 0.98 ± 0.02 are 
obtained in this case. 

In estimating the efficiency for the inclusive mode, relative contributions from all 
exclusive decay modes are fixed in the Monte Carlo simulation with the exclusive 
branching ratios obtained later in this analysis, giving an efficiency value of (64.05 ± 
0.03)%. The branching ratio for the inclusive mode is therefore obtained 

B{r~ ^ K-Xur) = (1.52 ± 0.04,^^%, (7) 
where X represents any system of neutral particles (including Kg^s). 

6.2 Selection for the K^-reduced and K^-enriched samples 

The analysis of exclusive modes involving kaons divides the one-prong r decays into K^- 
reduced and i^£-enriched samples, by using cuts on the two variables 50 and 6e. In the 
ir£-enriched samples, the net energy deposited in the ECAL should be at most a factor of 
two larger than in the HCAL, to remove the background hemispheres with unreconstructed 
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Figure 4: Fitted distribution for the inclusive one-prong r decay sample (top). The 
dots with error bars correspond to data. The fit, the expected pion ( including the small 
background from electrons and muons) and kaon contributions are indicated by the solid, 
dotted and dashed curves, respectively. The goodness- of-fit for the overall x-,,^ distribution 
is x^/^c?/ = 63/57. The residuals between the data and the fit are also shown in the 
bottom plot. 
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Figure 5: The distribution of charged kaon momenta in the inclusive one-prong r decay. 
The dots with error bars correspond to the fitted data and the histogram is for the Monte 
Carlo simulation. Errors are statistical only. 
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Mode 


fK{%) 


Nk 


X^/ndf 




fi(10-3) 


K- 


5.75 ±0.19 


2032 ± 66 


71.9/57 


11 ±4 


6.96 ±0.25 




1.92 ±0.10 


923 ± 47 


68.3/57 


15 ±5 


4.44 ±0.26 




1.55 ±0.23 


131 ± 19 


52.0/47 


3± 1 


0.56 ±0.20 




3.43 ± 1.09 


22 ±7 


33.5/35 


< 1 


0.37 ±0.21 


K-K^ 


9.20 ±0.96 


150 ± 14 


27.2/27 


2± 1 


1.62 ±0.21 




9.23 ± 1.18 


78± 10 


19.7/20 


1 ± 1 


1.43 ±0.25 




8.89 ±4.44 


4±2 


3.2/6 


< 1 


< 0.18 (95% C.L.) 



Table 2: Summary of statistics, overall charged kaon fractions in the Xj^ fit, the test of 
goodness- of-fit for the overall x^^ distribution in each mode, estimated qq backgrounds, and 
the branching ratios obtained later. Errors are statistical only. 



7r°'s. Meanwhile, the requirements of 6(j) < 0, 6e > 1.0 and 50^ ± > 12 are placed on 
the events to reduce the non-Kl background. Also, K~ K^n^^n^) candidates should be 
consistent with r decay, i.e., M{K~7r^) in K~ K^ti^Vt or M(A'~7r°7r°) in K~ K^tt^tt^i't- less 
than [Mt- — M^o). To further suppress the charged pion background from the 7r~7r°(7r°) 
modes, the angle ttopen between the direction pointing from the interaction point to the 
energy- weighted barycentre in the HCAL and the 7r'^(7r'^) momentum vector, is required 
to be smaller than 6° in selecting the K~ K^n^{n^) candidates. Finally, the K~ K^n^ir^ 
candidates are required to have Mlir^n^) greater than 0.5 GeV/c^, to reduce the K~K^7r^ 
contamination due to the decay 7r°7r°. 

The above selection criteria define the i^^-enriched samples for K~K^, K~ K^tt^ and 
K~ K^TT^Tc^ in the h~ , h~TT^ and h~7T^7T^ modes. The ii'^-reduced samples, involving decay 
modes like K~ , K~t^^ and K'tt^tt^, are defined as those which do not satisfy these criteria. 
Finally, the kaon fraction for each sample is determined from the fit to the corresponding 
Xtt distribution. 

6.3 Charged kaon fractions in the iC^-reduced samples 

The ii'^-reduced samples , h^ir^ and H^tx^tx^ are used to measure the branching ratios 
for K'Ut-, K^tt'^Uj. and K^tx^tx^Vt-. Fits to the distributions (Fig. ^ are performed in 
each momentum bin, in order to reduce uncertainties from dynamics in the simulation. 
The Xt^{tx) parameters in each momentum bin are directly taken from the values obtained 
in the fits to the inclusive mode as given in Table ^ Tiny fractions of muon and electron 
backgrounds are directly estimated from Monte Carlo and fixed in the fits. Table ^ gives 
the results of the fits. Figure |^ shows the charged kaon momentum distributions for 
both data and Monte Carlo simulations using the branching ratios obtained later. Good 
agreement is observed. 

For the H'tx^tx^tx^ mode, no cut is applied on the calorimetric measurement, and the 
kaon momentum is required to be greater than 5 GeV/c, reducing the pion background. 
The fit is only done for the whole momentum region because of low statistics. A generator 
with a. V — A matrix element and r phase space factor is used to produce Monte Carlo 
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Figure 6: The distributions in the K^-reduced samples. Data are shown with error 

bars. The fit results are shown as solid curves; the kaon and pion components are given 
by the dashed and dotted curves, respectively; the small electron and muon contributions 
are included in the pion component. 
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Figure 7: Kaon momentum distributions for the K^-reduced samples: (a) h~ mode; (h) 
/i""7r° mode; (c) /;,~7r°7r° mode. Data and Monte Carlo simulations are shown in dots with 
errors and in histograms, respectively. Errors on the measured kaon yields are statistical 
only. 
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Figure 8: Fitted x,r distributions in the K^-enriched samples: (a) for h^K^, (b) for 
h'K^ir^ and (a) for h~K^2TT^. The dots with error bars correspond to data. The fit 
function and the pion (including the small electron and muon background) and kaon 
contributions are shown in solid, dotted and dashed, respectively. 

events for the decay r~ — * K'TT^Tr^ir^Ur, giving an estimate of the average K/tt separation. 
A fit to the corresponding distribution is performed, and the result is given in Table |^. 
Figure P(d) shows the x-„ distributions and the corresponding fits in the i^~7r°7r°7r° mode. 

6.4 Charged kaon fractions in the A'^-enriched samples 

Decay modes such as K~K^, K~ K^t^^ and K~ K^tt^ti^ are measured by fitting the 
distributions in the i^^-enriched samples. 

In the analysis of h~K^ mode, the fits to the distributions in the i^^-enriched sample 
are also performed in momentum bins. This sample contains the tracks going through the 
ECAL cracks in order to enhance the selection efficiency, consequently introducing some 
electron background which is taken into account in the fits and has no influence on the 
kaon content. The fitting procedure is the same as used above, yielding the results shown 
in Table ^ and Fig. p. The charged kaon momentum distribution in the h~K^ mode is 
also shown in Fig. ^, indicating a good agreement between data and the Monte Carlo 
simulation. 

Fits to the full momentum range are applied for the K^tx^ and h^K^TT^ir^ modes 
because of low statistics. To reduce the pion background, an additional 5 GeV/c cut is 
employed for the minimum kaon momentum. The average K/n separation is obtained 
from the simulations. Table ^ and Fig. § show the results of the fits. 
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Figure 9: Charged kaon momentum distributions for data (dots with error bars) and the 
expected signal from simulation (histogram). 

6.5 Branching ratios for decays with charged kaons 

The measurement of the kaon yields for all exclusive modes is used to determine the 
corresponding branching ratios after applying the efficiency matrix obtained from Monte 
Carlo (see Table |^). In the simulation, form factors for the two- and three-meson decay 
channels are derived from Ref . [ 0, |l^ , and the four- meson r decay channels, for which 



no model is available at present, are generated according to the general V — A weak matrix 
element folded with phase space distributions for hadrons. 

All channels studied in this analysis are correlated through vr" reconstruction and 
contamination. Their branching ratios should thus be computed globally, by solving 
the seven linear equations with the efficiency matrix given in Table ^, in order to take 
the correlations properly into account. Since the channel K~t]i>t- can also feed into 
each studied mode, the corresponding subtraction is performed using the measurement 
B{t- K-riu^) = (2.9iij) x lO^^ from ALEPH [0. The final branching ratios for 
one-prong r decays involving charged kaons are given in Table 

7 Channels involving K^^ only 

This section studies the K'^n^ , K'^tt^tt'^ and K'^tt' decay modes, in which the 

potentially large backgrounds expected from the decay channels 7r~z/^, n'TT^Ur and 
n^n^n^Ur are separated by means of a fit method based on the previously defined variables 
6e and 6(p. 



14 



Mode 


R- 


R-n'-> 




R-TT^n'^TT^ 


R-R'' 


R-R"7t" 


R-R^^TT^ 


R-ri 


R- 


65.37 


4.25 


0.34 


- 


15.24 


1.28 


0.14 


0.28 




0.42 


44.26 


7.41 


1.18 


4.69 


8.89 


2.86 


1.84 


R-TT^TT^ 


~ 


0.50 


23.32 


12.40 


2.90 


3.69 


6.30 


5.34 


R-nW* 


- 


- 


0.28 


8.10 


- 


1.26 


2.82 


1.82 


R-R'^ 


0.55 


0.66 


- 


- 


16.86 


1.82 


0.48 


0.26 


R-R'>7r'> 


- 


0.48 


0.63 


- 


1.32 


10.04 


3.74 


0.20 


R-R\K^' 


- 


- 


- 


0.42 


- 


0.62 


4.54 


0.24 



Table 3: Efficiency matrix (in percent) for the exclusive r decay modes involving charged 
kaons. The K^-reduced samples involve the first four decay modes, while the K^-enriched 
samples include the last three. The generated decay modes are given in the first row, and 
the reconstructed modes in the first column. 



7.1 Selection for the samples 

In general, non-i^'^ background can mimic a signal through (i) h~ with upward 
fluctuation of calorimetric energy; (ii) 7r° leakage into the HCAL and (iii) unreconstructed 
7r°'s in the ECAL. As done in selecting the i^'^-enriched samples for the purpose of 
measuring the charged kaons, the selection criteria are based on the measurements from 
both the HCAL and the ECAL. The candidates are required to have at least 5 GeV 
energy deposited in the HCAL, to remove the background from a high energy 7r° leakage. 
The net energy deposited in the ECAL must be smaller than 15 GeV and less than two 
times the energy deposited in the HCAL, rejecting background with unreconstructed vr'^'s. 
According to the Monte Carlo study, the 7r~7r° or tt^tt^tt^ backgrounds may produce a fake 
signal in the HCAL when part of vr" energy leaks into the HCAL through the cracks 
between ECAL modules. To remove this background a quantity (pmoduie, defined as the 
azimuthal angle for the energy-weighted barycentre in the HCAL, is required to be outside 
a 5° window centred on each ECAL crack. A further suppression of this background is 
achieved by requiring that 6e be larger than 1 for K^tt~ and 1.5 for K^tt~tt^{tt'^), at the 
cost of a few percent signal loss. The charged hadron momentum is required to be less 
than 30 GeV/c for K^n-, 25_GeV/c for_^7r-7rO and 20 GeV/c for KOtt-tt^tt^. To be 
consistent with r decay, the K'^tt^tt^ or K'^ir'n^n'^ candidates must satisfy M(7r~7r°) or 
M(7r~7r°7r°) less than {M^- — Mj^o). The angle Oopen defined in Section 6.2 is required to 
be smaller than 8° for the K^ti^ti^ mode and 5° for the K^ti^ti^ti^ mode. The above cuts 
select three samples for the study of the K^7r~, K'^n^n^ and fC°7r~7r°7r° modes. 

7.2 Measured fractions 

The two variables 6e and 6(f) are used to distinguish the events from the background. 
In previous analyses [ H, H, the Kf^ signals are selected with cuts on 6e and 6(f), while 
in this analysis a fit method is applied. It is observed that systematic effects from the 
dynamics are reduced, and at the same time the distributions of 6e and 6(f> in data and in 
the simulation can be compared, providing a way to estimate the systematic uncertainties. 
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Rather than performing a two-dimensional fit, a one-dimensional fit is exploited using the 
variable ^ 

e = -^('5i?-A + (50) (for 5^, > 0, 50 > 0), (8) 

and 

e=^(5i?-A-50) (for 5s > 0, 50 < 0), (9) 

where A represents the 6e cut used for selecting the sample. This procedure is 
equivalent to projecting each event onto one of the two ^ axes shown in Fig. |iy together 
with scatter plots for data and Monte Carlo. Large 50 and 6e offsets are observed for 
the K^n~{K^) events, while they are relatively smaller in the K'K^Itc^) sample, because 
of the harder charged kaon momentum. The h~{7r^) backgrounds contribute to the 
enriched sample because of an upward energy fluctuation, yielding a long tail in the 6e 
distribution. The definition of C, is optimized for the K^Tr~{K^) signal. Events with a 
positive (negative) sign of 50 always correspond to a negative (positive) sign of ^. As a 
result, the majority of K^^s are distributed in the region of large positive and are thus 
efficiently separated from the non-i^T^ background. 

To measure the K^^ fractions, the shapes of ^ distributions for both the modes 
and the non-ZT^ modes are parametrized from the Monte Carlo simulation, using a set 
of Gaussian functions. Since the majority of K^^ events are on the positive ^ side, the 
negative ^ distribution can be used to check the calorimetric calibration and the simulation 
as well. 

In the fits to the ^ distributions, the means and the widths of all the Gaussian functions 
for the non-/r£ modes are left free in order to check the simulation and the calibration. 
The channels involving a K~K^ pair are also taken into account by incorporating the 
shapes from the simulation and fixing the fractions from the independent measurements 
of the branching ratios given in Table 0. The channels involving two neutral kaons cannot 
be separated from those involving a single neutral kaon because the two showers 
overlap in the calorimeters. A small difference is indeed observed in the simulated ^ 
shapes, but the fit is not sensitive to this effect. Table || gives the results obtained in 
the ^ fits. Figure |ll] shows all the ^ distributions and the corresponding fits to the three 
/T^-enriched samples. 



7.3 Branching ratios for decays with one only 

The relevant branching ratios are obtained from the measured yields using the 
efficiency matrix given in Table ||. Since the K'^ K'^ti' {Ti^)i'r channels can feed into the 
channels with a single neutral kaon, the subtraction of these backgrounds is necessary. 
The contribution of the Ki^K^n' mode is computed using the measurement of the 



KgKgTT'Ur channel, assuming CP invariance. The latest measurements [ H, |T3| of the 
KgKgiT' mode give an average (0.24 ± 0.05) x 10~^, significantly smaller than the value 
(0.75 ± 0.38) X 10~^ [ |14| used in the previous ALEPH analyses [ |], §], which relied on the 



assumption B{t~ — > K^K^^tv~ Vr) / B{t~ K^K^tx^Vt-) = 1/4. After correcting for this 
effect, the published measurements [|^, ^ are in good agreement with the value obtained 
in the present analysis. Since the K'^K^TT- mode can also give a contribution to all the 
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Figure 10: The scatter plots (5e vs 64>) in h mode for (a) data, (h) Monte Carlo 
events for t~ — > K^t:~{K^)ut- (r^5x statistics of data), (c) Monte Carlo events for 
T- K-K^{'K^)v^ (r^ 8x data), and (d) h-{'K^)v^ 5x data). The two ^ 

axes considered in the analysis are drawn with straight lines as shown in each plot. 
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Figure 11: The ^ distributions for the sample (a) in h~ mode; (h) in h~Tr^ mode; 
(c) in /i~7r°7r° mode. Data are given with error bars. For ^ > 20, the bin width is 

increased to 2.0. The fits, signals (h~ = 7r~ ) and r backgrounds are shown as solid, dashed 
and dotted curves, respectively. Even though the distributions for positive and negative 
are disconnected and vanish at ^ = by construction apart from binning effects, it is 
convenient to display the overall distribution in the same histogram. 
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Mode 








^ (non-K9) 


-0.04 ±0.03 


0.00 ± 0.03 


-0.03 ± 0.16 


(non-Er£) 


0.98 ±0.02 


0.96 ± 0.02 


1.00 ±0.08 


fKl (%) 


15.24 ±0.68 


8.67 ±0.76 


8.12 ±3.41 


Nki 


937 ± 41 


299 ± 26 


17 ±7 




72.5/69 


46.0/41 


14.0/21 




2±1 


16 ±6 


6±2 




25 ±4 


12 ±3 


< 1 


B (10-^) 


9.28 ±0.45 


3.47 ±0.53 


< 0.66 (95%C.L.) 



Table 4: Results from the fits to ^ distributions. The ^ parameters for the non-K^ 
modes, the fraction of K^, the corresponding number of 's, the test of goodness-of- 
fit in each mode, the estimated backgrounds from K^K^n~7c^ and qq, and the branching 
ratios obtained later. 



Mode 


Kill" 












45.03 


4.20 


0.23 


32.11 


3.86 




2.34 


24.42 


1.17 


8.71 


2.93 




0.24 


1.16 


8.64 


0.53 


3.46 



Table 5: Efficiency matrix (in percent) for the exclusive t decay modes involving one 
neutral kaon only. The feedthrough backgrounds from the channels involving two neutral 
kaons are taken into account, as detailed in the text. The generated decay modes are given 
in the first row, and the reconstructed mode in the first column. 



final states as shown in Table ^, the subtraction is performed using the branching ratio 
B{t- K^Kln-u^) = (1.01 ± 0.26) x 10"=^ given by ALEPH [ |]. As for the decay 
mode K^K'^TT~7c^, a first measurement of the branching ratio of (0.31 ± 0.12) x 10^^ 
for the KgKln~TT'^ part has been obtained by ALEPH [ 0. For this mode, the ratios 
KgKg : K^K^ : K^K^ of 1/4 : 1/4 : 1/2 are assumed, yielding the relevant contribution 
to each final state given in Table ^. Finally, the branching ratios are computed using the 
efficiencies given in Table ^ and the numbers of neutral kaons, after doing the K^K'^t^'t:^ 
and qq background subtraction, where the branching ratios for are given by assuming 
equal and contributions. The results are given in Table ^. 

8 Systematic uncertainties 

The sources of systematic uncertainties in this analysis are divided into several parts: 
general selection criteria, dE/dx response, calorimetric measurements, feedthrough 
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Source 


Sel 


dE/dx 


Laiori 


Bkg 


MC stat 


Dynam 


Total 


K X 


U.o 


2.0 




0.7 


0.4 


U 


2.7 


A 


n A 
U.4 


i.o 




u.y 


U.D 


n 
U 


o n 
z.U 




0.7 


4.8 




2.0 


1.2 





5.4 




4.1 


12.1 




23.5 


2.3 


5.8 


27.5 




15.8 


4.1 




22.1 


4.8 


10.0 


29.6 


K-K^ 


0.6 


1.4 


3.2 


1.8 


2.7 


5.0 


6.9 




0.9 


10.7 


4.2 


6.9 


4.2 


4.0 


14.6 




0.6 




2.5 


2.2 


1.4 





3.7 




1.1 




5.4 


6.1 


3.4 


6.0 


10.8 



Table 6: Summary of systematic uncertainties (relative in percent). See text for details. 



background, Monte Carlo statistics and decay dynamics. The estimates are given in 
Table ^ and are detailed below. 

8.1 Selection efficiency 

In this analysis, one-prong r decays are classified into h' , h^n^, h^-K^-K^ and h^n^Ti^TT^ 
samples. The uncertainties on the selection efficiencies, studied in Ref. [ H, correspond 
to the general r-pair selection, the lepton/hadron separation, the 7/7r° recognition, the 
handling of fake photons and the hadronic interactions. In the general r-pair selection, a 
0.2% uncertainty is given for all non-K° channels, while a 0.3% uncertainty is assigned to 
the channels involving K^. An uncertainty of 0.2% on the lepton/hadron separation 
is estimated for all relevant channels, i.e., without involving vr^'s. The systematic 
uncertainty related to vr*^ reconstruction increases with the number of vr^'s involved. Errors 
of 0.6%, 2.6% and 15.7% are estimated for the channels involving one, two and three 
7r°'s, respectively. Concerning fake photons, a detailed investigation was presented in 
Ref. [ |]. It was found that there are some discrepancies between data and Monte Carlo, 
resulting in corrections which are included in the efficiency matrix in Table ^ and Table |^: 
(-1.49 ±0.32)% for K'Ur, (-0.58 ± 0.40)% for R-n^r, (+4.28 ± 3.17)% for /sf-TrVz/,, 
(+0.21 ± 1.49)% for ir-vrWz/^, (-1.62 ± 0.54)% for h'K^Ur, (+1.06 ± 0.93)% for 
K'^tt't^^i't- and (—0.59 ± 0.61)% for K~ K^tt^Ut-. The errors are taken as the systematic 
uncertainties. All the above four contributions to the general selection criteria are put 
together and given in Table |^. 

8.2 dE/dx measurement 

Three systematic effects are relevant to the dE/dx measurement. Following Ref. the 
first part is the pion calibration. The corresponding uncertainty is estimated from the 
study of the inclusive mode, in which the fits to the distributions give two values for 
the error on the kaon fractions (see Table |I]), the first obtained from the three-parameter 
fits by leaving the x,r(7r) parameters free and the second from the one-parameter fit when 
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fixing the X7r(vr) parameters. The pion calibration uncertainty is obtained from their 
quadratic difference, giving a 0.07% uncertainty on the absolute kaon fraction. This value 
is applied in a relative way to all channels involving charged kaons, yielding 2.3% {K~X), 
1.4% {K-), 4.3% (K-ttO), 11.1% (i^-TrV), 3.2% (K-ttOttOttO), 0.9% {K-R'') and 10.2% 
{K~ K^ir^), respectively. 

The second part is related to the shape of the pion resolution function F^(x^) 
derived from the muon samples. The largest contribution to the uncertainty comes 
from the statistical errors in the parametrization of F^(a;,r); which correspond to 
relative uncertainties of: 0.7% (K'X), 0.4% (R-), 1.2% (ir"7r°), 3.2% (K-7r°7r°), 0.9% 
(ir-7r°7r°7r°), 0.3% {R-R^) and 2.9% {R-R^-k^), respectively. Another contribution 
arises from the energy dependence of the parametrization using r~ — >■ ^'u^^Ut and 
77 ^~ events. The corresponding uncertainties are: 0.1% {R~X), 0.1% {R~)-, 

0.8% (K~vr°), 1.7% (K-7r°7r°) and 0.3% (R^R^), respectively. The other modes with 
charged kaons are selected with a 5 GeV/c minimum momentum cut, thus reducing the 
effect studied above to a negligible level. 

The third part comes from the R/ vr separation, which follows the calibrated velocity 
dependence obtained from pions. There is no additional calibration error for kaons in 
principle. However, for convenience the kaon resolution function, Eq. (^), is constructed 
from Monte Carlo tracks, hence generating small statistical uncertainties of 1.0% {R~X), 
0.8% {R-), 1.4% (R-TT^), 3.4% (ir-vrV), 2.4% (ir-vrOvrV), 1.0% {R-R'') and 1.4% 
{R~ R^ir^), respectively. These errors are translated into the systematic uncertainties for 
each channel, and finally added to the other two dE/dx uncertainties in quadrature as 
shown in Table |^. 

8.3 Calorimetric measurements 

The systematic errors due to the uncertainties in the calorimetric measurements are 
investigated by comparing the 6(f) and the 6e distributions in data and Monte Carlo 
simulation. From the statistical errors in the calibration procedure discussed in 
Section 5.3, an uncertainty is derived, through the 6(f) cuts, on the efficiencies for the R~R^ 
and R'R^n' modes, and the feedthrough background from other r decays, yielding 
systematic uncertainties of 1.0% and 1.1% the R~R^ and R~R^ti~ modes, respectively. 
Negligible uncertainties affect the other modes. 

The precision of the 6e measurement gives corresponding uncertainties of 2.9% 
(R-R"^), 3.9% (ir-ifV), 0.2% {R-) and 0.2% {R-n''). 

In selecting the i^^^-enriched sample, a cut on the ratio of energy deposited in the 
ECAL and the HCAL is used to remove the background from /i~-7r° and /i~7r°7r°. To 
investigate the possible biases, two r decay samples h~ and Sh" are studied. In the h~ 
mode, events are required to have a positive 6(f> and a minimum 5 GeV HCAL energy. 
The comparison between data and Monte Carlo gives a (—2.7 ± 1.0)% correction to the 
efficiency estimated from the Monte Carlo simulation. In the Sh" mode, a similar study is 
performed, yielding a (— 0.2±0.7)% correction. In this analysis, two hadrons are produced 
(charged hadron + R^) and an average of these two results yields (—1.5 ± 1.0)% for the 
R^h~{Tx^) modes, leading to a correction included in the efficiency matrix in Table ^. 

In obtaining the ^ shape for the R^ signal, an uncertainty may arise from the validity 
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of the calibration constants obtained from the non-i^T^ samples when applying to the 
sample with overlapping showers. To evaluate this effect, the 3/i r decay mode is used to 
mimic the signal: 5e and 50 are computed, using only one charged track momentum 
but the full calorimetric measurement for the three charged hadrons. As a result, the ^ 
distribution looks very similar to the signal, providing a relevant comparison between 
data and Monte Carlo simulation. Applying the derived correction to the ^ fit yields a bias 
^fxl — (0.06 ± 0.16)% for the fraction in the K^t:~ mode. A relative uncertainty on 
the corresponding branching ratio is therefore estimated to be 1.6%. The same correction 
is also applied to the K'^n~n'^ mode, giving a 3.9% relative uncertainty on the branching 
ratio. 

Finally, the systematic effect due to the calorimetric calibration is investigated. 
Analogous to the systematic uncertainty study for dE/dx measurement, this effect is 
studied by fitting the mean and the resolution of ^ {i.e., three-parameter vs one-parameter 
fits), giving 1.6% and 3.4% relative uncertainties on the branching ratios for K^n'Vr and 
K^Tr~7r^ i/t, respectively. 

8.4 Background subtraction 

Since all one-prong r decay with kaons are studied simultaneously, except for K~r] and 
the modes with two neutral kaons, the uncertainties due to the feed-across effects are 
included in the branching ratios, taking into account the errors from the fits and the 
errors on the branching ratios used. 

Background from p and Oi can feed into the K'^7r~ and K'^n'n^ modes when the tt^ 
leaks through an ECAL crack. The cut on (pmoduie is designed to remove this background 
(see Section 7.1). However, the Monte Carlo simulation produces a somewhat broader 
shower shape in the calorimeters, leading to an uncertainty in calculating the barycentre 
in HCAL. To investigate this effect, a special sample is selected by requiring (pmoduie to 
be in a 5° window centred on an ECAL crack and ^ > 5, enriching the background from 
p and fli up to 50% of the total. The comparison between data and Monte Carlo shows 
that this background is over-estimated by a factor of (1.2 ±0.1) in the signal region. This 
correction is taken into account in obtaining the ^ shapes used in Section 7.2 for the /i~7r° 
and /i~7r°7r° modes with unreconstructed 7r°'s. The corresponding errors are treated as 
systematic uncertainties, yielding 0.7% for K'^tt'i't and 1.0% for K'^7r~7r^UT-. 

Turning to the non-r background, mainly qq, the uncertainty arises from the dynamics 
to produce kaons. To investigate this effect, tighter cuts are applied on the hemisphere 
opposite to the one studied, rejecting hemispheres with a hadronic invariant mass greater 
than 3.0 GeV/c^. More than 50% of the qq background is removed at a cost of about 
1.2% in efficiency. Applying these cuts to the data reduces the fitted charged kaon yield 
by 85 ± 15. The Monte Carlo predicts that the kaons from qq background arc reduced 
by 40 ± 6, while those from r decay decrease by 47 ± 3, resulting in a loss of of 87 ± 7 
kaons which agrees well with data and indicates that kaon production in these very low 
multiplicity qq events is well simulated within a 40% uncertainty. The same uncertainty 
is also given to the channels with neutral kaons. 
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8.5 Monte Carlo statistics and dynamics 

Statistical uncertainties in the determination of the efficiency matrices hsted in Table ^ 
and Table |^ are taken into account. However, the most important source of uncertainty 
originates from the decay models used in the generator, affecting both the kaon and the 7r° 
momentum spectra. This effect may produce biases when determining the momentum- 
dependent K/t^ separation, the 7r° reconstruction efficiency and the identification 
efficiency. All of these contributions are discussed below. 

For the K~ , K'n^ and K^7r~ modes, the decay dynamics is well understood, and the 
corresponding uncertainties are neglected. Since the inclusive K~X mode is dominated 
by the well established K~ and K~tt^ modes, and the ffis to the distributions are 
performed in momentum bins, the uncertainty due to the unknown dynamics is neglected. 
In the K~K^ and K~ti^tt^ modes, the ffi to the x-,^ distribution is performed in momentum 
bins, reducing the uncertainty on the determination of K/n separation to a negligible level. 

In order to estimate the effect of and 7r° momentum spectra, the efficiencies 
are studied as a function of the hadronic mass. A mass uncertainty of 100 MeV/c^ is 
assigned for the modelling of the K~K^, K'K^n^, K^ti^-k^ and K^ti'ti^ decay modes. 
The variations in the efficiencies are taken as systematic uncertainties, yielding 5.0%, 
4.0%, 5.8% and 6.0% for the above decay modes, respectively. For the K'tt^h^tt^ mode, 
a 10% uncertainty is assigned for the unknown dynamics. 



9 Investigation of mass spectra 

Mass spectra are investigated in the K'n^, K^tc^tt^, K^ti^ti^tt^ , K^n~ , K^7r~n^ and 
K^K~T[^ modes. Because of the overlap between the hadronic showers, the K~K^ mode 
cannot be studied. The overlapping problem also affects the investigation of the total 
invariant mass for K'^'k~tt^ and K^K~tt^ . Nevertheless, the study of subsystems such as 
7r~7r° and K~tt^ can still provide some useful information about the relevant dynamics. 
A more complete assessment of the resonance structure taking into account all isospin- 



related channels is available f 15 . 



9.1 Mass spectra in the decays involving one K only 

The study of mass spectra for the K~7c^, K^tt^ti^ and K'tt^tt^tt^ decay modes is different 
from the conventional method, in which cuts are applied to isolate the charged kaon at a 
cost in efficiency. Instead, this analysis generates the mass plots by fitting the number of 
charged kaons from the distribution in each mass slice. 

Since the minimum kaon momentum in the decay K^n^Ur is about 3.5 GeV/c, 

a cut is set for excluding the low momentum tracks in the following analysis. A one- 
parameter fit is performed to extract the number of charged kaons in each mass slice. 



giving the mass plot shown in Fig. |T2|, with a 14% background level from other r decays 
involving charged kaons and by construction no background from pions. Except for a 
small excess on the high mass tail in the data which is discussed elsewhere [ the 
K*{892)~ signal is well predicted by the Monte Carlo simulation, with the background 
mainly distributed on the low mass side. 
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Figure 12: The K^tt^ invariant mass spectrum for the decay channel K'tt^Ut-: a) 

linear scale, b) log scale. Data are shown as dots with error bars and the Monte Carlo 
predictions are shown in the open histogram together with the background distribution in 
shaded. 
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In the K~iT^TT^ analysis, the requirement on the minimum kaon momentum is raised 
to 5 GeV/c. The mass spectra for the ii'~7r°7r°, K'n^, and vr^vr'' systems are investigated 
accordingly. Figure |T3| shows these mass spectra for the K^tt^tt^ mode. Many r decay 
channels feed into this final state and contribute about 60% of the total. However, since 
most backgrounds are distributed in the region of low mass, some useful information 
related to the channel r~ — ^ K'lr^ir^i'T- can still be deduced. A peak at the i^*(892)~ 
mass confirms the expected K*{892)^tt^ intermediate state. In general, i^'i(1270) and 
-ft'i(1400) can be produced in the KtttiVt- channels. The model in Ref. [ |TU[ assumes 
-ft'i(1400) dominance. However, the K~Tt^TT^ mass spectrum in Fig. ^ is consistent with 
an admixture of these two states. 

Because of low statistics in the i^~7r'^7r°7r° mode, the investigation of the mass spectra 
is difficult. However the vr^'Tr^Tr'^ invariant mass can offer some useful information regarding 
the importance of the K~ri channel. Only 5 ± 4 out of 22 ± 7 events are found in the Svr 
mass region (0.4-0.65 GeV/c^), ruling out a dominant r] contribution, in agreement with 
measurements in the 77 — 77 mode [ O. Most of the charged kaons are found in the 



region M(37r'^) above 0.65 GeV/c^, in which signal and background from K K^n^Vj. with 
Kg —> 7r°7r*^ overlap. 

9.2 Mass spectra in the decays involving one only 

Candidates for the channels involving one are isolated by using cuts on both 50 and 
^. In the K^n~ mode, the charged pion momentum can be small, resulting in a gap 
between the and vr" showers. The requirements 6(f) < —10 and ^ > 10 remove almost 
all non-K^ background, yielding a purity of 86% for the decay channel r~ K^T^'Vr. 
The flight direction is defined from the colliding point to the barycentre of the energy 
deposited by the in the HCAL. The distribution of the invariant mass M{K^tt~) is 
shown in Fig. |14|, with a clear i^'*(892)~ signal seen. The backgrounds from non-A'* decay 
modes, such as the channels K^K^t[~ Vr and K^-K^-K^Vri give a 10% contamination, while 
the qq background contributes 4%. The Monte Carlo simulation is observed to be in 
good agreement with data, showing that the energy measured in the HCAL is adequately 
simulated. 

In the K^Ti~Ti^ mode, the total invariant mass provides less information because of the 
rather poorly reconstructed momentum and the low statistics. However, the study of 
the TT'TT^ invariant mass distribution can give useful constraints on the relevant dynamics. 

The K^TT~TT^ candidates are selected by requiring the total measured energy in the 
HCAL to be greater than 10 GeV, the angle aopen (defined in Section 6.2) less than 6° 
and ^ > 5. Under these additional requirements, the backgrounds from p~ /a^, i^*(892)~, 
K'^K^T[~ and K^K^ti^ti^ are suppressed to levels of 14%, 26%, 8% and 8%, respectively. 



The TT TT invariant mass is shown in Fig. |15|. An incoherent sum of a p Breit-Wigner 



signal and the shape of the K vr refiections is used to fit to the background-subtracted 
'7r~7r° invariant mass plot, giving (72 ± 12 ± 10)% for the Kp fraction, where the second 
10% uncertainty is due to the background subtraction. 
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Figure 13: The invariant mass spectra for the K^tt^tt^ mode in data (dots with error bars), 
Monte Carlo predictions (histograms) and the expected background distribution (shaded 
histograms). The model predictions assuming iiri(1400) dominance [ are shown in 
dashed histograms. 
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Figure 14: The K^t:~ invariant mass spectrum for the decay channel t~ K^-K'Ur- Data 
are shown as dots with error bars. Monte Carlo predictions and the expected background 
distribution are shown in the open histogram and in the shaded one, respectively. 




Figure 15: The 7r~7r invariant mass spectrum in the decay channel t~ K^tt'tt Ur- 
Data are shown as dots with error bars. The Monte Carlo prediction (open histogram) , the 
background contribution from p~/ai (hatched histogram) and other t backgrounds (shaded 
histogram) are shown, while the signal follows the dashed histogram. 



27 




Figure 16: The K^tt^ invariant mass spectrum for the K^K~tt^ candidates. Data (dots 
with error bars), Monte Carlo prediction (open histogram), background contribution from 
K^7i^{7i^) (hatched histogram) and other t backgrounds (shaded histogram) are shown, 
while the signal follows the dashed histogram. 



9.3 Mass spectra in the decays involving a K pair 

As mentioned before, the total invariant mass distributions in K~K^j^Vt- and K~ K^j^t^^Vj. 
cannot be usefully studied. However, the K~7i^ mass plot in K~ K^h^Ut can be 
investigated. A requirement < —2 is added for the K^h~Ti^ candidates, leading to 
a purity of 65%. The K~n^ mass plot is shown in Fig. 16 with a clear K*{892)~ signal 
observed. Although the feedthrough K^n^i^n^) background can fake a K*{892)' signal, 
it cannot account for such an excess in the K*{892)^ mass region. In the model of Ref. [ 
10| used in the Monte Carlo simulation, the K*K^ intermediate state in K~ K^ti^Vj. is 



treated as strongly suppressed, which disagrees with this observation. After subtraction 
of the background, the K~7t^ invariant mass distribution for the decay r~ — > K^K^n^Ur 
is fitted with two components, K*K and pvr, showing that the K*K intermediate state 
contribution in K^K~Tx^Vr is dominant (larger than 86%). This observation is consistent 
with the result from the K~^K~tt~ mode studied in Ref. [ |l|. 



10 Results and discussion 

All measured branching ratios in this analysis are listed in Table ^. The inclusive 
branching ratio for r~ —>■ K~Xvt^ is measured to be (1.52 ± 0.04 ± 0.04)%, where X 
represents any system of neutral particles. The sum of all measured exclusive branching 
ratios gives (1.46 ± 0.05 ± 0.04)%, in agreement with the inclusive value. 

The (i^7r)~ mode is studied in two final states, K~tt^ and K^-k' , both showing 
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Decay 


B (10-3) 


T- K-Xur 


15.20 ±0.40 ±0.41 




6.96 ±0.25 ±0.14 




4.44 ±0.26 ±0.24 




0.56 ±0.20 ±0.15 


K^T^^T^^T^^Vr (excl. rj) 


0.37 ±0.21 ±0.11 




1.62 ±0.21 ±0.11 




1.43 ±0.25 ±0.21 




< 0.18 (95% C.L.) 




9.28 ±0.45 ±0.34 




3.47 ±0.53 ±0.37 




< 0.66 (95% C.L.) 



Table 7: Summary of branching ratios obtained in this analysis. For the inclusive mode, 
X represents any system of neutral particles. 



dominance of i^'*(892) . 

The measurement of r decay into K^tx'^ti'^u^ completes the study of the (i^vrvr)" final 
states started in Ref. [ |, |]. The branching ratio B{t- K-n^n^Ur) = (0.56 ± 0.20 ± 
0.15) X 10~3 is significantly smaller than the value (1.1-1.4) x 10^'^ predicted in the model 
of Ref. [ |ll|, but in good agreement with the value 0.4 x 10"^ in Ref. [ |16]. This channel 



suffers from a large background contamination as evidenced by the study of invariant mass 
spectra. Despite this, a dominant K*{892)^tt^ intermediate state is observed. Background 
from K' found to contribute to the K tt^tt*^ mode, resulting in a difference 

from the previous ALEPH measurement B{t- R-ti^ti^p^) = (0.80±0.20±0.20) x 10"^ [ 
wherein this background was ignored. 

Assuming that the decay K^-k^-k^Vt proceeds via incoherent intermediate states 

Kp and K vr, one can derive the following contributions: 

B{t- (K\)-Ur ^ K^n-n^Ur) = (0.97 ± 0.44 ± 0.36) x 10'^, (10) 
B{t- ^ (Kp)-Ur ~^K^n-7i^Ur) = (2.50 ± 0.57 ± 0.44) X 10^1 (11) 

These results are in agreement with the measurement based on the Kg mode [^. 

Finally, it is found that the decay r~ K~ K^tt^v^- mainly proceeds via the K*K 
intermediate state. 



11 Conclusion 

One-prong r decays involving either charged and/or neutral kaons are measured in this 
analysis. As already done in the other two ALEPH analyses [ [1|, 0], the study of final 
states is extended up to four hadrons involving kaons, giving a complete measurement for 
all sectors of r decays involving kaons. The results are given in Table |^ and also shown 
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in Fig. [T^ to compare with other experiments [ [T^ . Agreement is observed with all the 
published data. A forthcoming paper will combine all ALEPH measurements on r decays 
involving kaons and will give the relevant physics implications [ |l5l . 
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Figure 17: The published branching ratios for one-prong r decays with kaons The 
black dots correspond to this analysis. The detected modes are given in parentheses 
for each experiment. 
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